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Abstract 



We construct local M-operators for an integrable discrete time version of the classical Heisenberg 
magnet by convolution of the twisted quantum trigonometric 4x4 R- matrix with certain vectors in its 
"quantum" space. Components of the vectors are identified with r-functions of the model. Hirota's 
bilinear formalism is extensively used. The construction generalizes the known representation of M- 
operators in continuous time models in terms of Lax operators and classical r-matrix. 

> 

^ ' 1 Introduction 



The unified treatment H of non-linear soliton equations as hamiltonian systems having enough number of 
conserved quantities in involution is based on the (classical) r-matrix. Its role is to provide a universal form 
of Poisson brackets for elements of Lax operators. An alternative though less popular point of view on the 
r-matrix (which we are going to follow) comes from the Zakharov-Shabat approach El, (2) which consists 
in representing soliton equations as 2D zero curvature conditions for a pair of matrix functions (called L 
and M operators) depending on a spectral parameter. The r-matrix works there as a machine producing 
M-operators from L-operators. 

In the paper || , for the simplest example of the lattice sine-Gordon (SG) model, we have found a similar 
machine in the completely discrete set-up, i.e., in lattice integrable models with discrete time. Remarkably 
enough, it appears to be a quantum R-matrix with the "quantum" parameter q related to the time lattice 
spacing. The formula for local Af-operators in the discrete case has as simple structure as in the continuous 
one (see ( |l.5| ) below) with the i?-matrix in place of the r-matrix. This comes as the result of a computation. 
At the moment we do not try to " explain" why a typical quantum i?-matrix takes part in a purely classical 
problem. 

In this paper we show that the construction still works for a much more general model - partially 
anisotropic lattice Heisenberg magnet (HM) in discrete time. (This model contains discrete versions of 
SG, KdV and other equations as special cases.) Although the final result looks like a straightforward 
generalization of the one for the lattice SG model, the derivation gets much more involved and some important 
modifications are necessary. 

Let us recall the r-matrix construction of M-operators for continuous flows. Let Ci(z) be a classical 2x2 
L-operator on ID lattice with the periodic boundary condition Ci+n{z) = Ci(z); z is the spectral parameter. 
For the moment we assume the ultralocality - Poisson brackets between elements of the Ci{z) at different 
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sites are zero. The monodromy matrix %{z) is ordered product of L-operators along the lattice from the 
site I to I + N - 1: 

%{z) = Ci+n-^z) . . . £ l+1 (z)Ci(z) . (1.1) 

Hamiltonians of commuting flows are obtained by expanding logT(z) in z, where T(z) — tvTi(z) does not 
depend on I due to the periodic boundary condition. All these flows admit a zero curvature representation 
with the generating function of M-operators given by || , jjj 

M l (z;w)=T-\w)tr l [r{z/w)(Ti(w) <g> I)] . (1.2) 

Here r{z) is a (classical) 4x4 r- matrix, acting in the tensor product of two 2-dimensional spaces, tri means 



trace in the first space, I is the unity matrix. Expanding the r.h.s. of eq. (1.2) in w, one gets M-operators 
depending on the spectral parameter z. From the hamiltonian point of view, the zero curvature condition 
follows from Poisson brackets for elements of the L-operator. A similar r-matrix formula exists for M- 
operators in non-ultralocal models, though, in this case the r-matrix is not skew-symmetric. In general, 
Mi(z;w) is non-local. 

It is well known Q , |l[] how to construct local ALoperators from the generating function. Suppose there 
exists zo such that det£;(zo) = for any I. This means that Ci{zq) is a 1-dimensional projector: 

rl s Hi)){W)\ 

A(zo) = J ^ j (1-3) 

where 

<p\ = (a, a) (i-4) 

are two-component vectors and P(l) is a scalar normalization factor. Components of the vectors as well as 
the P(l) depend on dynamical variables. It is easy to see that Mi(z; zq) is a local quantity: 



(0(1) 


r(z/z Q ) 


a(l-l)) 


(M) 


a(l 


-1)) 



(note that the normalization factor cancels) . The scalar product is taken in the first space only, so the result 
is a 2x2 matrix with the spectral parameter z. It obeys the zero curvature condition 

d t Ci{z) = Mi +l (z)Li{z) - Ci{z)Mi{z) . (1.6) 
Let us outline the results of this work. We consider the completely discretized partially anisotropic 



(XXZ) HM and represent M and L operators Mi(z), Ci(z) by formulas of the type (1.5) 



M i( z ) = / R n\\ n — ' (I- 7 

(/3(l)\a(l - 1)) 

{p{l)\a(l - 1)) 

In the r.h.s., R(z;q,£), R(z;q,^') are quantum 4x4 i?-matrices (to be specified below) with the "quantum" 
parameter q and the Drinfeld twist parameters related to the space-time lattice. The vectors |a(Z)), 
are ^ e same as in eq. ( [1.5| ), 



It is worth noting here that the continuous and discrete time models have a common L-operator; formula 



(1.8) gives its i?-matrix representation. 

To visualize these formulas, it is convenient to use the language of the algebraic Bethe ansatz [|| , 0] . 
The scalar product is taken in the "quantum" (vertical) space, so one gets a 2x2 matrix in the "auxiliary" 
(horizontal) space: 
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(P\R(z)\a) 



The M-operator (1.7) generates shifts of a discrete time variable m. The discrete zero curvature condition 

Ml+l, m (z)£l, m (z) = Cl, m+ l[z)Ml,m(z) (1.9) 

gives rise to the HM equations of motion in discrete space-time. 

The key ingredient of the construction is to pass from the initial dynamical variables to the pair of vectors 
|ck(0)i |/3(0)- I n the papers || on exact lattice regularization of integrable models components of these 
vectors were expresed in terms of canonical variables of the model. Those formulas looked quite complicated 
and were hardly considered as something instructive. Here we shed some light on their meaning. Using 
equations of motion of the completely discretized model (derived from the discrete zero curvature condition 



(1.9)), we show that (suitably normalized) components of the vectors |a(l)), |/3(Z)) are T-junctions. 

The r-function is one of the most fundamental objects of the theory (see e.g. ||, ||). It is worth stressing 
that quantum i?-matrices in the classical discrete problems can hardly be noticed until one reformulates the 
non-linear equations and elements of the L-M pair entirely in terms the r-functions. That is why we are 
to make a long excursion into Hirota's bilinear formalism jl(|-[jl2| (Sect. 2). Our guiding principle is Miwa's 
interpretation |13|| of the discrete time flows, in which discretized integrable equations are treated as members 
of the same infinite hierarchy as their continuous counterparts. A general method to produce discrete soliton 
equations was developed in juj, where, in particular, the discrete isotropic (XXX) HM model was proposed. 
Following these ideas, we deal with bilinear form of the discrete XX Z HM model. The matrix L-M-pair for 
the latter is derived in Sect. 3 from basic linear problems for a scalar wave function |Q, juj, |fL5[ , 

Like in the SG model, there are two lattice versions of the classical HM: the discrete HM model on a 
space-time lattice |l^| and the model on a space lattice with continuous time introduced in |fl6f for the more 
general XYZ-ca.se. It is interesting to note that both of them were proposed back in 1982 but the ideas 
underlying one and another seemed to be " orthogonal" and did not intersect until the very recent timeQ. A 
motivation of this work was just an attempt to understand one of the models in terms of the other one. 

The manner of exposition here is different from the one in S, where we started with either Faddeev- 
Volkov or Izergin-Korepin L-operator for the lattice SG model in terms of lattice fields and passed to the 
r-function by means of a special substitution. In this paper we start directly from the r-function (and bilinear 
equations it obeys) rather than from lattice spin variables. The reason is twofold. First, the equations of 
motion in terms of lattice spin variables are too complicated and, anyway, are of no practical use for us. 
Second, in the discrete case there is no canonical way to introduce spin variables whereas Hirota's bilinear 
formalism is gauge-invariant and free of artifacts. Having this in mind, let us remark that by Heisenberg 
model we mean here nothing more than a convenient name for the properly reduced 2D Toda lattice hierarchy 
p2| . (Actually we deal with a model which is slightly more general than the XX Z HM itself.) In this sense 
a name like partially anisotropic chiral field model would be also appropriate. Needless to say that in the 
discrete set-up the specific features of each of the two models become irrelevant || . What is relevant is the 
type of reduction of the 2D Toda hierarchy that is the same in both cases. 

Treated as a hierarchy, the bilinear equations imply infinitely many discrete variables (commuting flows) 
from the very beginning. Generally speaking, any two of them could be chosen as space-time coordinates. 
However, for our purpose we need a distinguished space-time lattice. In Sect. 4 we introduce such a lattice 
with coordinates denoted by I, m and derive the matrix L-M-pair for translations I — > Z + l, m — > m+1. These 

1 On the quantum level, these ideas were partially linked together in [p^|-[po|, where quantum SG model in discrete space-time 
(a quantization of Hirota's discrete SG equation) was constructed; a generalization to afBne Toda field theories on the lattice 
was suggested in 
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L and M operators are represented in the form (1.8), (1.7) in Sect. 5. In Sect. 6 we discuss the continuous time 
limit of the construction and show how eq. ( |1.5| ) is reproduced. As it was already mentioned, the r-matrix 
in that formula turns out to be not necessarily skew-symmetric. This fact is a signal of non-ultralocality of 
the model in general case. Sect. 7 is a conclusion where we point out a few unsolved problems motivated by 
our results. 



2 Hirota's bilinear formalism 



In this section we present the results of the papers [13J, |14j] in the form convenient for our purpose. We 
illustrate formulas by the graphical representation of discrete flows suggested in the review [23 to which we 
refer for more details. 



2.1 General form of 3-term bilinear equations 

The key object of Hirota's approach is r-function. The r-function r = r(a, 6, c, . . .) (as a function of the 
discrete variables a, 6, c, . . .) obeys a number of bilinear partial difference equations. To each discrete variable 
a a complex parameter A a £ C ("Miwa's variable") is associated. One may think of the A a as a continuous 
"number" which marks the discrete flow. For any two discrete flows a, b we put 

A a fc = A Q — Af, , A a b = —A(, a ■ 

Let us recall how to compose bilinear equations. Each triplet {a, b, c} of discrete flows gives rise to a 
3-term bilinear equation [^2[, [fL3| for the t: 



Ab c r(a + 1, b, c)r(a, b + 1, c + 1) 
+ A ca r(a,&+ l,c)r(a+ l,b,c+ 1) 
+ X ab r(a, b, c + l)r(o + 1, b + 1, c) = . 



(2.1) 



All other variables which the r may depend on enter this equation as parameters. Each quadruplet {a, 6, c, d} 
gives rise to another 3-term bilinear equation: 



A Q d-Ab c T(a + 1, b, c, d + l)r(a, b + 1, c + 1, d) 
AbdA ca T(a, b + 1, c, d+ l)r(a + 1, b, c + 1, d) 
A c dAabT(a, 6, c + 1, d + l)r(a + 1, 6 + 1, c, d) = . 



(2.2) 



Remark Links between eqs. (2.1), ( |2.2| ) exist in both directions. On the one hand, eq. (2.1) is a particular 
case of (2.2) when Ad — * oo. (According to [13], this limit means that the dependence on d in the r disap- 
pears.) On the other ha nd, e q. (pT^), though linearly independent of eqs. (|2.l|), is an algebraic consequense 
of equations of the type (2.1) written for the triplets {a, b, c}, {a, 6, d} and {a,c, d}. In this sense all what 
we are going to derive in the sequel follows already from eqs. (2.1). 



2.2 Basic bilinear equations for the discrete HM model 

Let us consider five discrete flows and denote the corresponding variables by p, p' , n, u, v. We denote 
the r-function by t(j>,p' ,n,u,v) = t^ p (u,v). In the latter notation the u,v are separated from the others 
because they will play the role of chiral space-time coordinates. We call them chiral variables. 
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Bilinear equations for the r-function of the XX Z HM model are obtained from (| 
the reduction 



p+l,p' + l ; 



(U,V) =7^f 1 ('U,u), 



(2.2) by imposing 
(2.3) 



so we are left with four independent variables. Let us choose them to be p, n, u, v. 

Remark The discrete KdV equation is a particular case A„ = oo, so the dependence on n disappears. The 
discrete SG model corresponds to further specification A P ' = A p . 

Using the graphical representation introduced in |23) , we can display the above configuration of discrete 
flows in the figure (the graph of flows): 

A p ' 




The dashed line is drown here to indicate that the reduction (2.3) is the same as in the case of the ID Toda 



chain in discrete time (a detailed discussion of this point see in |23|]). Indeed, the vector field dp = —d n + d p 
defines the flow p corresponding to the dashed line; in these terms the reduction acquires the more familiar 
Toda-like form t p+1, p +1 = t p, p. This turns into a differential condition as A p — > oo, \ p > — > A n . In this case 
one gets the isotropic (XXX) HM model in discrete time. 

Remark Looking at the graph makes it easier to deal with different discrete flows when there are many of 
them. We refer to |^3| for the exact meaning of the graphical elements. Here we only note that keeping in 
mind solutions of the finite-gap type (see e.g. 0), (24j|), the reader may think of this figure as drown on a 
patch of a Riemann surface with local coordinate A . Miwa's variables are coordinates of punctures on the 
surface. The lines are then cuts between the punctures, which give rise to discrete commuting flows on the 
Jacobian via the Abel mapn. 



The condition (2^3) allows one to g et rid of p' if necessary. Let us give a list of bilinear equations which 
are obtained in this way from eqs. (2T), ( |2.2|) . In front of each equation we indicate the triplet or quadruplet 
which it comes from. 



The simplest equation involves the variables p, n only: 

{pp'n} : \ p , n T%~l(u,v)T%±[(u,v) - \ pn rP +1 (u^)^ 1 (u,v) = \ p>p (t%(u, v)) 2 . 



(2.4) 



It is the discrete time ID Toda chain equation in bilinear form [|12|. In our case it plays the role of a 
constraint on the dynamical variables since the space-time coordinates u, v enter as parameters. 

Equations of the next group contain shifts of p, n, u only: 



{pp'u} : 



r p+i 



(u, v)t* + 1 (u + 1,v)- \ pu tI + \ (u, v)tP +1 (u + 1, v) 



= Vp t £( u + MK +1 (u,v). 



(2.5) 



{pnu} : A u „t£ +1 (u,v)t£ +1 (u + 1,v) + X pu t^ +1 (u,v)tP +1 (u + l,v) 

= X pn r^u+l,v)^+l(u,v), (2.6) 



For the general theory of finite-gap solutions to Hirota's difference equation see 
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{p'nu} : \ un TP(u,v)TP +1 {u + 1,v) + \ p , u tP +1 (u,v)tP(u + l,v) 



(2.7) 



{pp'nu} : \ p „\p> u tP +1 (u,v)t% 1 (u + 1, v) - X P > n X pu T^ + l(u, v)t^ 1 (u + 1, v) 

= \p>p\ u nT%{u,v)T%(u + l,v) . 

Similar equations can be written for p,n,v - it is enough to replace u by v everywhere. 



Remark Combining eqs. (2.5), (2.6), 



eq.^ 



cq-(M 



\ p > u tI + \{u+1,v) X vn T% +1 (u + l,v) 
and plugging ( p7z| ) into the l.h.s., we obtain eq. ( |2.8| ). 
Equations involving both coordinates u, v read: 



(2.8) 



{nuv} : \ un T%(u,v + l)i% +1 (u + l,v)-X vn T%(u + 1, v)t% +1 (u, v + 1) 

= \ uv t%(u + 1,v + 1)t% +1 (u, v), 



(2.9) 



{pp'uv} : \ p > u \ pv tP +1 (u,v + 1)t% +1 (u + 1, v) -A pu A p ^T^ +1 (u, v + 1)t£ +1 (u +1,v) 



(2.10) 



{puv} : A p „< +1 (u, v + 1)tP(u + 1,v)- X puT P +1 (u + 1, v)t%(u, v + 1) 

= \ uv tP{u + \,v + 1)tP +1 {u,v), 



(2.11) 



{p'nuv} : -\ p , u \ vn TP +1 (u,v + 1)tP(u + l,v) + \ p , v \ un TP +1 (u + l,v) T P(u,v + 1) 

= X p ^XuvT^I(u + 1,v + 1) T p +1 (u,v), (2.12) 



{p'uv} : \ p , vT p +1 (u,v + 1)tP +1 (u + 1,v)-\ p , u tP +1 (u + 1, u)t£ +1 (u, v + 1) 

= X uv tP +1 (u + 1,v + 1)tP +1 (u,v). 



(2.13) 



{pnuv} : -A p „A lJn r^ +1 (u,'W+ 1)t£ +1 (u + l,v)+X pv X un T% +1 (u + 1, v)tP +1 (u, v + 1) 

= VLrffu+M + l)^^!;), (2.14) 

The above list of linearly independent 3-term bilinear equations is by no means complete. The full list 
contains many other equations which either follow from the already written ones or can be derived from 



"higher" analogues of eqs. (2.1), (2.2) for more than 4 variables (which, in their turn, are algebraic corollaries 
of eqs. (2.1)). Some of them are given in the next subsection. 
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2.3 Some useful corollaries of the basic equations 



Note that the pair of equations ( fHl| ), (respectively (pl^) , can be considered as a linear 

system for tP +1 (u, v + l)r^(u+ l,v), t^ +1 (u + 1, v)t^(u, v + 1) (respectively, for rf l+1 (u, v + 1)tP +1 (u+ 1,v), 
Tn+i( u + ^-i v ) T n +1 { u i v + !))• Resolving these systems, we get 



\'vKnT^{u + 1,V+ 1)tP +1 (u,v) + \ pu \'nT^t\{u + 1,V + 1)t^ +1 (u,v) 

AtP +1 (u,v + 1)tP(u + 1,v), 



(2.15) 



A p '„A to tP(m + l,v+ 1)tP +1 (u,v) + \ pv Xp'nT^t\{u + l,v+ 1)t% +1 (u,v) 
At%(u,v + 1)t% +1 (u+1,v), 



(2.16) 



\v^unT^ +1 {u + 1,V + 1)t^ +1 (u,v) + \p> u ^pnT%(u + 1,V+ 1)t^+ x (ti, «) 

ArP(u,v + l)rP +1 (u + l,v), 



(2.17) 



where 



Ap„A TO r^ +1 (u + 1, v + 1)t^ +x (w,«) + Ap'„A p „T^(u + l,i> + l)r^ x (u,u) 
A7^- 1 («,t; + l)7* +1 («+l J i;) J 



(2.18) 
(2.19) 



These equations are equally useful in what follows. They are e ven more in formative than eqs. ( 2.11 )-( 2.14 ) 
since allow one to "fuse" the variables, i.e. to put X u — X v . (In (2.11 )-(2.14) this leads to the identity = 0.) 
In this case t{u + 1, v + 1) converts into t(u + 2) and we get: 



A„„A p ' u tP +1 (w - 1,v)t*(u + 1,v) + A pu Ap'„r^ +1 (u - l,u)r^ x (u + l,v) 

(Xp'nXpn ~ X 2 un )TP{u,v)TP +1 {u,v) , 



(2.20) 



Kn\ u T? l+1 (u - 1,v)t% +1 (u + 1,v) + \pnXp> u tV +1 {u - 1,v)t%(u + l,v) 

= (Xp'nXpn - Xl n )TP +1 (u,v)TP +1 ( U ,v) . 

The same equations can be written down for v in place of u. 



(2.21) 



At l ast, we present two other useful equations which arc obtained from ( |2.5| )-( |2.8| ), ( 2.15 )-( 2.18 ) and 



( 2.20 ), ( 2.2 lh (together with their w-counterparts) by a procedure similar to the one explained in the remark 
after eq. (^.§|): 



XpnXp'nT^^U + 1,V + 1)T% +1 (U,V) - X un X vn T%(u + 1, V + 1)t%(u,v) 



(2.22) 
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\ p > u \p< v \pnT% 1 (u + l,v + 1)tP +1 (u,v) - XpuXpvK'n^l-iiu + l,v+ 1)t% + {(u,v) 
= \ p , p At%(u + 1,vK(u,v+1). (2.23) 

To conclude the section, let us mention that the complete list of linearly independent equations for the 
r-function (even restricted by the 3-term ones) is somewhat longer. Here we have selected those used in the 
sequel. 

3 Linearization of the discrete HM model 

The bilinear equations from the previous section can be represented as compatibility conditions for an 
overdetermined system of linear problems for a "wave function" ^ . This is what we mean by the linearization. 
For us this is a systematic way to find L-M-pairs. 



3.1 Scalar linear problems 



The bilinear equations ( |2.4|) -( p7l4"l) follow from compatibility of a system of linear equations for a "wave 
function" * = *£> p ' (u,v). The prototype of linear equations for * is Q, Jl|], |1| (see also j|| for a 
review): 

*(a + 1,6) = *(a, b + 1) - A ab ) ; \ \ w _ '- ■ ,{ ^(a, b) 



r(a + l,b)r(a,b + 1) 
where a, b stand for any two elementary discrete variables. 



The reduction condition (2.3) for the \l/-function reads 
where z is a spectral parameter. Hiding the variable p' with the help of this condition, we have: 



{nu}: *E(«+l) = tt£ +1 («)-A 



r p n {u)T p n+l (u+l) 
''tI +1 (u)tI{u + 1) 



{pu} : ^ n +1 (u) = + 1) - V^StO^*"^) > 



(3.1) 

(3.2) 

(3.3) 
(3.4) 
(3.5) 



(and similar equations for v in place of u). In (3.3)-( p.5| ), v is supposed to be the same everywhere and, 
therefore, skipped. These equations are basic tools for deriving matrix L-M-pairs. 

3.2 Vector linear problem 

Combining equations (3.3)-(3.£), one can represent translation of the vector I ^p+i/ u \ J along the u- 

. (3.6) 



direction in the matrix form 
/ **(«+!) \ / 



7 2 \ 

6 /\pn 



( U ) T n+l( U ) 



pn p+ i 
Th, 



2" + Ap' u — — 



^(uJtT^u+i) y 



<s 



Let us change the gauge passing to the wave function 



$i(u) 
$ 2 (w) 

where 2?(u) is the diagonal matrix 

Then the linear problem ([T^) acquires the form 

*i(«+l) 

$2(W + 1) 



V(u) 



\lnziTP(u) \ 

xJz^rP+^u) ) 



L {u \z) 



*i(«) 
$ 2 (w) 



The L-operator L^ u \z) can be compactly written in terms of the three fields 



T n+l(«> W ) ' 



ip (u, u) 



T% +1 (ll, v) 



For brevity we also use the notation <fi(u,v) = \ip°(u, vj] 2 . The L-operator reads 



L {u \z) 



("+!) 

0(«) 



i>-(u + l) \ 



(u)4>(u + 1) 



Note that 



z\U{u)^u + l)^{u) X p , u -^)+z*^±^ 

4>{u + l) <p(u) J 



detL^ (z) — Xp U Xp> u — X un z 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



We call the L^- W \z) chiral L-operator because it shifts the chiral variable of the wave function via ( |3.8| ). In 
the next subsection we study the discrete zero curvature condition with matrices of the type ( [3.1C ) . 



3.3 The zero curvature condition for chiral L-operators 

Let 



■ C=(u,v+1) 



D=(u+l,v+l) 



A=(u,v) 



B=(u+l,v) 



be an elementary cell of the u, ii-lattice^]. In this notation (borrowed from p7[-[p!9[) the L-operator ( 3.10| ) 
reads: 



'In general the coordinate axes are not orthogonal to each other; in particular, at X u = A„ the lattice collapses to a ID one. 
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-(«) 

J B^A 



(z) 



x pu 



HA) 



' 4>{A)^B) 



(3.12) 



Similarly, we introduce another chiral L-operator, L^]_ A (z), which is given by the same formula with 4>(C) 
in place of </>(-£?) and X pv , X p > v in place of X pu , X p ' u , respectively. 

The discrete zero curvature condition 



L 



d^b( z )Lb^a 



is equivalent to the following non- linear equations of motion for the fields ip , "0 



o „/,±. 



(3.13) 



^m {C)^\D)[^ + {C)^-{C) -tP + (B)^(B)} , 



(3.14) 



X pv i> {D)^ (C)^-(B) + X P , U ^{A)^(C)^-{D) 
X pu ^(D)^(B)ir{C) + X p/v ^°(A)^ a (B)^-(D) , 



(3.15) 



A jro V°(^)V' (C)V + (C) + Vu^°(>l)^ (CO^ + (A) 
X pu iP°(D)ip (B)rJj+(B) + A p ^° W° 



(3.16) 



Equations of motion for the discrete HM model of XXX type given in 1 14 are written for a different choice 
of dynamical variables. It would be useful to establish a direct correspondence between them. 



Remark The discrete KdV equation in the Faddeev-Volkov form |17[, faM is reproduced from ( 3.15 ) or ( 3. 16 ) 
in the limit A„ — > oo when the n-dependence disappears, so that tp~ = ip + — 1 and we are left with the 
single field The Faddeev-Volkov chiral L-operators for the discrete KdV equation are reproduced from 



(3.12) in the same limit provided the renormalizad spectral parameter £ = zX pn is finite. 



3.4 Antichiral L-operators 

Here we introduce another type of L-operators which will be called antichiral. 



Let us put An = A n in the chiral L-operator (3.12) and, correspondingly, identify u with n. In this way 
we get the operator which generates the translation 



A=(u,v,n) — > A' = (u, v, n + 1) 



in the 3D lattice spanned by u, v, n: 



Xpn 



H^) 



HA) 

zx} n <t>(A)<t>(A^+(A) 



iX 



i Q4t) \ 
V(A)0(AT) 



Xp'n 



cf>(A^) + <j>[A) ) 
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V )0V i 



p+i 



»pn p 



i 'n+l 
A P" p+1 



(^) 



+ A 



r p n+ MKtl(A) J 



(3.17) 



The similar notation (£p = (u + 1, v, n + 1), = (u + 1, u, n — 1) etc) will be used for the vertices in the 
upper and lower layers of the 3D lattice. 



Now we define the antichiral L-operator to be 



L 



Apn Ap'n 



L 



(v) 



(z). 



It generates a discrete flow v (which we call antichiral) defined by the vector field <9 5 
of variables. We can display the flow v on the graph of flows from Sect. 2.2: 



A,; 



A (1 



A u 







• 


\ v 


n 





(3.18) 

-d n + d v in the space 



Computing the r.h.s. of ( [3.18 ), we get: 



L ci^a( z ) 



( 



\ 



A # ^ 



He 1 ) 



<f>(Cl) " pn " pv 0(A) 



zXpn 



A A 



(3.19) 



(^-counterparts of eqs. (2.6), (2.7) have been used). 

Remark Although we do not need the antichiral flows themselves, we found it instructive to mention them 
here in order to motivate the choice of the L-operator in the next section. Let us stress that antichiral flows 
have "equal rights" with the chiral ones. The symmetry can be seen from the figure and from the explicit 
form of the L-operators. Note also that in the case of the discrete KdV or SG models there is no difference 
between chiral and antichiral flows. 



4 "Composite" L and M operators 

The chiral and antichiral L-operators introduced in the previous section are building blocks for more com- 
plicated ones obtained as their ordered products. Not all of them have an i?-matrix representation of the 
desired form. Our next task is, therefore, to find a special pair of flows l,m such that the corresponding L 
and M operators meet the requirement. The lattice spanned by l,m will be called the space-time lattice. It 
is embedded into the 3D lattice with coordinates u, v,n. 

Our experience in the lattice SG model suggests that the good candidates are "composite" L and M 
operators which generate translations along diagonals of the chiral space-time lattice. In the discrete HM 
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model this prescription works literally for the M-operator only. That is why, contrary to the tradition, we 
discuss the M-operator first. The X-operator requires an important modification coming from the fact that 
in general the space-time lattice is embedded into the 3D lattice in a different way than the chiral one. 



4.1 The "composite" Ai-operator 

Consider the "composite" operator which transfers the wave function along the diagonal C — > B: 



Mb^c(z) = z L {X vn z z - X pv X p * v )L)£_ A (z) L£!_ A (z) 



From (3.12) we have: 



-(«) 



(4.1) 



M B ^ c (z). 



■J>(B) 



<t>{C) 



i jr(D) \ 



(C; B^) zX vn Jjy^ Z Xp' u Xp v 



<Am(C, B) = Xp' u ^ | - X p >.. 



4>{C) 

<KB) 



(4.2) 



The diagonal elements are brought into this form using eqs. ( |2.6| ) and (2.7). The non-diagonal elements are 
obtained with the help of eqs. ( |3.15| ), ( 3.16 ). 



Let us express the matrix elements in terms of the r-function. For any two vertices X, Y of the lattice 
we set 

W{X,Y)^[r^\X)rl +1 {X)r^\Y)rl +1 {Y)Y • 
Using eqs. (pl|) , (|t|), we find: 



M 



B^C 



(*) = 



1 



/ zX un rl +1 {C)rl +1 {B) 
-z-% u X p , v tI +1 {B)tI +1 {C) 



W{C,B) 



_\ \ 2 _P+1 

^uv^pn < n _|_ 



\{A)t%+\D) 



-X uv X^ n rP +1 (A)rP(D) 



\ 



zX vn Tl +1 {B)rl +1 {C) 
-z-% Iu X pv tI+\C)tV +1 {B)) 



(4.3) 



Moreover, with the help of eqs. ( 2.17 ), ( 2.1 8| ) it is possible to express the matrix elements through the 
r-functions at the points A and D only, while the dependence on the points C, B is confined in the prefactor 
VF" 1 (C, B). Indeed, the r.h.s. of eqs. ( 2.17 ), (|2.18| ) are just the two terms on the diagonal of the matrix in 
(4.3). Substituting them by the left hand sides, we see that the matrix depends on A and D only. This is 
used in Sect. 5. 



4.2 The "composite" L-operator - first version 

The first candidate is suggested by our experience in the discrete SG model 

£dJ-a( z ) = z 1 ^'dI-b( Z )^ j B^a( Z ) ■ 



(4.4) 



Using the bilinear equations from Sect. 2 many times, we can express its matrix elements through the r- 
function: 



1 



XuvW(A,D) 



( A uv ) 

r (uv) 
\ L "2X 



r (uv) 
*-12 



C 



(uv) 



(4.5) 
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with 



47 } = \Mz)tI +1 (B) T v + \C) - \ pu Uz)T p n+1 {C) T v + \B) , 



^ = X? n z (Cu{z)rlX\{B)r^\C) - Uz)< + + \{C)r^ {B)) , 



(4.6) 
(4.7) 
(4.8) 



+ T^T" U P ' v XunUzKtl(B)rP(C) - \ p , u \ vn C v (zKXl(CM(B)) . (4.9) 

/Ann /\-)ir» ' 



Here 



(4.10) 



As it will be clear later, the desired i?-matrix representation for this operator does not exist since the 
matrix elements have a "wrong" dependence on z. We should seek for an L-operator which would depend 



on z as in (4.3) 



4.3 The "composite" L-operator — improved version 

The above L-operator can be "improved" by an additional translation in n. The right choice is the operator 
which generates the translation 

A = (u, v, n) — ► = (u + 1, v + 1, n — 1) 
in the 3D lattice spanned by u, v, n. Its elementary cell is shown in the figure: 

C = (u, v + 1, n) D = (u + 1, v + 1, n) 



A = (u, v, n) 



B = (u + 1, v, n) 



C l = (u,v + l,n- 1) 



D L = (u+ 1, v + l,n - 1) 



A* = (u, u, n — 1) 



L? 1 = (w + l,t>,n- 1) 



One of possible ways to represent such an L-operator in terms of already mentioned ones is: 

i -l 



(4.11) 
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with L^l_ nl (z) a s in Jo]) and C^ Az) as in ( j4.5| ). Matrix elements of /^^^(z) are Laurent polynomials 
in z. From (4.11), ( 3.1 2| ) and (4.6)-(4.10) we see that the diagonal elements are of at most 3-d degree, the 



right upper (resp., left lower) element is of at most 4-th (resp., 2-nd) degree. 

It turns out that this L-operator actually has a surprizingly simple structure - much simpler than one 
could expect from (4.11). In particular, using the bilinear equations, one can show that the highest degrees 
of z cancel leaving us with at most z and z^-terms in the diagonal elements and z-independent terms in 
the non-diagonal ones. 



This cancellation could be expected from the following equivalent representation of the L-operator (4.11) 
which makes use of the antichiral flow introduced in Sect. 3.4: 

C m ^ A (z) = L { ^\^ ci {z)L { ^l^ A {z) . 



This means that the flow generated by the L-operator ( 4.11 ) is the "superposition" of the chiral and antichiral 
ones. 

Computing the L-operator by means of any one of the two formulas, we find: 



(z\ X HA) 



C D i^ A (z)- 



>(L>i) 



(4.12) 



>(L»i) 



HA) ■ 



In terms of the r-function we have: 



/ zX un X vn TP(D)TP +1 (A) 

+ 



W(A,Dl) 



+ Z 1 \puXp V \p' n Tn1il(D)TP +1 (A) 



AAp n r^ +1 (B)t% +1 (C) 



z\ pn \ p , n rl + _\{D)rl +1 {A) 
-z-\ n X p , u X p , v TP(D)TP +1 (A) J 



(4.13) 



(A is given in (2.19)) 



4.4 The space-time lattice 

Let us consider the 2D lattice spanned by coordinates I, m which are flows generated by the L and M 
introduced above. The vector helds 



di = d u + d v - d n 



f) = f) _ a 



are then "coordinate axes" of this lattice. (Note the equivalent representations of these vector fields in terms 
of antichiral flows: di = d u + d v = d u + d v , d m = d u — <9 S .) Shifting the origin of coordinates, we always can 
embed this lattice into the 3D one as the 2D plane defined by the homogeneous linear equation 



u + v + 2n = 



(4.14) 



provided (it, v,n) £ Z 3 . In the sequel this 2D lattice is refered to as the space-time lattice (STL). Note that 
the chiral lattice is the plane n = 0, i.e. it is the linear projection of the STL along the n-direction. 



14 



The space-time coordinates I, m are introduced by the formulas 



I = -(u + v), 



m = —(u — v) 
2 V ' 



(4.15) 



provided the point (u, v, n) belongs to the STL. (Under the same condition I could be equivalently defined 
as I = —7i.) 

Remark After an obvious matching of the notation, the zero curvature condition on the elementary cell 



of the STL for the constructed "composite" L-M-pair takes the form (1.9). It gives rise to a system of 
non-linear equations for fields at the vertices of eight neighbouring cubs of the 3D lattice. This system is a 
direct consequence of the basic equations of motion (3.14)-( 3.16| ). 



5 i?-matrix representation of the L-M-pair 



Our goal is to represent the L and M operators (4.3), ( 4.13 ) constructed in the previous section as convolu- 
tions of a quantum R- matrix with some vectors in its "quantum" space. 



5.1 The quantum i?-matrix 

Consider the following quantum i?-matrix with the spectral parameter z: 



R(z) 



= i \2a(z) + b+(z) + b-(z)] I ® I + ^ [2a(z] - b+(z) - b-(z)] a 3 ® cr 3 
+ T [&+0) ~ b -( z )} (I® 0-3-03®!) + \c\z) {01 ® 01 +0-2® 0-2) 



( o(z) 



V 





b+(z) 
c(z) 






c(z) 





\ 



a(z) J 



(it,: are Pauli matrices), where 

a(z) = qz — q~ 1 z~ 1 , b±(z) = £ 1 (z — £ _1 ) , c(z) = q — q^ 1 . 



(5.1) 



(5.2) 



When necessary, we write R(z) = R(z;q,£); q is the "quantum" parameter and £ is the parameter of 
Drinfeld's twist |§, @ : 

R(z;q,Z)=F(Z)R(z;q,l)F(Z), 

The i?-matrix ( p.l\ ) satisfies the quantum Yang-Baxter equation for any q, £: 

Ri2{zi/z2)Ri3{zi/z z )R 2 3{z2/zz) = R23(z 2 / z 3 )R 13 (zx/ z 3 )R 12 (zi/ z 2 ) , (5.3) 

where Ri 2 (z) acts in C 2 ® C 2 ® C 2 as the R(z) in the first and the second spaces and as the / in the third 
one (similarly for i? 13 (z) and i? 2 3(z))- 



The first (resp., second) space of tensor products in ( |5.1| ) is called "quantum" (resp., "auxiliary") space. 
It is convenient to represent the R- matrix as a 2x2 block matrix in the "auxiliary" space. Let i,i' number 
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block rows and columns and let j,j' number rows and columns inside each block (i.e., in the "quantum" 
space). Then matrix elements of the i?-matrix (5.1) are denoted by R{z) 1 ^,. 

Let | a), |/3) be two vectors in the quantum space (see ( |l.4| )). Each block of the i?-matrix is an operator 
in the quantum space. Consider its action to \aj and subsequent scalar product with (j3\. The result is a 
2x2 matrix in the auxiliary space: 



(P\R{z)\a) w =Y^R{z%a r l3 j . 
jj' 



Substituting the matrix (|5.l|), we find: 



(p\R(z)\a) = 



Piaia(z) + f3 2 a 2 b + {z) 

( z{q[3 1 a 1 + £,(i 2 a 2 ) 
-z'Xq^Pxai + £/3 2 a 2 ) 



(3 2 aic(z) 
/3i«i6_(z) + (3 2 a 2 a(z) 



V 



(q-q 1 )Pia 2 



{q-q )@2ai 



-Z- 1 q- 1 r 1 {q^i+^2a 2 ) J 



(5.4) 



Hereafter all scalar products of the type (/3\R(z) |a) are understood as to be taken in the first ("quantum") 
space. 



5.2 The L-M-pair in terms of the f?-matrix 



We are going to bring the L and M operators (4.12), (4.2) into the form (5.4). The best result is achieved 
after the diagonal gauge transformation 



Mb^c(z) — > Mb^c(z) = 

t D i^ A {z) — > C D i^ A (z) = 
and shift of the spectral parameter 



tV + \B) T v +1 {B) 



M b ^c{z) , 



C D i^ A (z) 



z * kz . k — 



^pu^p' u 



(5.5) 
(5.6) 

(5.7) 



The matrices Ms^cikz), C D i i _ A (kz) are to be expressed through the r-function as is explained at the end 
of Sect. 4.1. 

In order to present the result, let us introduce the vectors 



k(u, v, n)) = 



\(3(u,v,n)) = 



HtP +1 (u,v) 



where 



^pu 
^pn ^p' u 



(5.8) 



(5.9) 
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and identify the parameters as follows: 

11 i 

(^vn^pu^p'u \ 2 j. / ^un^pu^p'v \ 2 j./ / ^un^vn^pu^pv \ 2 in\ 

T \ \ ) ' S — I T T \ ) i S — ~~ I x 2 \ \ ) (O.IU) 

^un^pv^p' v / \ ^Dn-^p^^p'u / \ p' u^p' v / 

(note that ££' = — /i 4 ). For brevity we also set 

P — [Apn Autt, Xvn^pU Apt; Ap'y, Ap'?;] 2 , 

then q — q~ x = — Ap n A u ^Ap _1 (recall that A is given by ( |2.19| )). 

After these preparations, the comparison with (5.4) immediately yields the formulas 

(p(B)\R(z;q,Z')\a(C)) = -Xuvp- 1 r^ +1 (A)^ +1 (A)C D i^ A (kz) , (5.11) 



(f3(D)\R(z;q,0\a(Al))=Ap- 1 T*+ 1 (C)iZ +1 (C)M B ^c(kz). (5.12) 



To recast them into the form (|L8|), ( |l.7| ) respectively, we, first of all, fix the constant time slice m = of 
the STL. It is embedded into the 3D lattice as the ID lattice with vertices at the points Ai = (I, I,— I), 
I 6 Z. Let us introduce the similar notation for other points: Bi = (I + 1,1,-1), Bi = (I + 1,1 — 1,-1), 
D\ = Ai +1 = (I + 1,1 + 1,-1 - 1) etc and put 

C l (z) = C Dl (z), M l (z)=M Sl ^ Al (z). (5.13) 

I 1 

Now we notice that 

< + x{u, v)t^ (u, v) = (VA^Ay)* {Ku + i, v , n) | a(tt _ i, w> n + i)) (5.14) 



due to eq. (2.20). This relation allows us to rewrite eqs. (5.11), (5.12) in the form 
where 

\m)) = = 1/3(5*)) 

V ^T*_f{i + i,i) ) V / 

and 7 is a constant: 

7 = [Aim Xpv Xp'v] 2 (Ap7^Ap' n ^un) ■ 

These formulas differ from the ones announced in the Introduction merely by the irrelevant prefactors. To 
match the notation, let us rename the vectors: |a(C/)) — > |a(0)i ~~ > 1/3(0); so horn now on we set 

/ ^{1,1 + 1) \ ( »T P _f{l + l,l) 

|«(0) = , |/3(0) = 1 • (5-17) 

V a*- 1 ^ 1 -(1,1+1) / V M-'^a+i.o 

Their location in the 3D lattice is shown in the figure: 
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The constant time slice m — is the chai n of vertices (. . . , Ai-i, Ai,Ai+±, . . .). The dashed line shows the 
translation generated by the M-operator ( 5.16 ). 

Remark The L-opcrator Ci(z) ( |5.15 ) has two degeneracy points: z = k and z = kq^ 1 . At the first one the 
i?-matrix is proportional to the permutation operator V in C 2 <g> C 2 : R(l; q, £') = (q — q^ 1 )!^. Therefore, 



i \a{l)) (j3{l)\ 

Ci{k) - A^„A , — — 

(lJ)T_ l+1 (l,l) 



(cf. (|l.3| )). At the second one we get: 



C l (kq- 1 ) = \$ n A- 



|/3(0) 



(5.18) 



(5.19) 



|/3(0) = 



o 

1 



|/3(0), 



«(0) = 



o 

mr 1 



Formulas of the type ( 5.15 ), ( 5.16 ) for the same L-M-pair can also be written in terms of the vectors |<5(0), 
1/3(0)- 



6 Continuous time limit 



In this section we show that the r-matrix formula for local M -operators (1.5) is a degenerate case of our 
R- matrix formula. It follows from Sect. 2 that X uv plays the role of lattice spacing for the discrete time 
variable m. At the first glance, the continuous time limit would then imply A„ — > A u , i.e. q — > 1, £ — > 1 
in the i?(z;g;£), so, in agreement with eq. (1.5), we do get a classical r-matrix. However, this would imply 



lim g _>i |/?;) = |az) that is certainly wrong in general. The correct limit to a continuous time coordinate is 



more subtle and needs some clarification. 



6.1 Limiting form of the M-operator 



We need a continuous time limit such that the L-operator (5.15) remains the same. The naive limit does not 
work because varying At, one changes the STL and, therefore, the L-operator itself. In the correct limiting 
procedure, the time lattice spacing should approach zero independently of \ u and At,. 
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To achieve the goal, we exploit the advantage of Miwa's approach [|13| and introduce w - another " copy" 
of the chiral flow v with Miwa's variable A^,. The bilinear equations from Sect. 2 apparently hold true for 
w and X w in place of v and A„ respectively. Now we can safely tend X w — > X u thus eventually getting a 
continuous flow, leaving A„ and all the other parameters unchanged. Set 



Xwn Xpu Xp f u 
Xun Xp W Xp' w 



and consider the limit e — > 0. We have 
9 = 1 



\p'u^un~^~ -^pu^p' n ^ ^ Ofs^) 



2 Xun ^pu A<p' u 



A-un Ap-u Ap' to 
Aiun Apiu Ap' u 



£ = i- 



Ap'xtA^n Ap U Ap' n 
2A<uri,Apii -t- Ap' u 



e + 0(e 2 ). 



(6.1) 



(6.2) 



Discrete M-operators are defined up to multiplication by a scalar function of z independent of dynamical 
variables. To pursue the continuous time limit, it is convenient to normalize the M-operators by the condition 
A4®(z) = I at £ = 0. Then the next term (of order e) yields local M-operator of a continuous time flow. 

To implement this project, consider the 4D lattice with coordinates (u,v,w,n). The m = slice of the 
former STL is embedded into this lattice as the ID sublattice with vertices Ai = (1,1,0,-1), I £ Z. The 
continuous time M-operator at the site Ai is obtained by expansion of the discrete M-operator Aig^^i ( z ) = 
X un k(z — z~ x )l + 0(e) which generates the translation Ai = (I, I, 0, —I) — > B' = (I + 1, I, — 1, —I). Passing 
to the normalized Af-operator, we have: 



M s ,^ A (z)=I + eM l (z) + O(e 2 ). 



(6.3) 



For clarity, the it, ^-section of the 4D lattice (v = const, n — const) is displayed in the figure: 




The point B[ tends to the point Ai as e — > 0, so the parallelogram collapses to the it-axis. Using eqs. (fO|), 
( |2.13 ), we obtain: 



Mi(kz) 



2Xp' U (z~ z ) 



1 , Is 

- - z+(z+-)Ui 

z z 



2[i 



^(1 + 1,1)^(1-1,1) \ 

T^\l,lV_ l+1 (l,l) 



2(1 



-2 T -l 



(1+1,1)^(1-1,1) 



v+1 



t\U)r P . l+1 (l,l) 



2^ + 2^ + z + --(z + -)U l 



where 



^\i + i,iV_ l+1 (i-i,i) 



u, = 



,(1,1)^(1,1) 



(6.4) 



Note the redundant freedom which allows one to add to the M(z) a term h(z)I with a scalar function h(z) 
independent of dynamical variables. This obviously does not affect the zero curvature condition ( |l.6| ). In 
the next subsection we shall use this freedom to redefine the M-operator. 
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6.2 Comparison with the r-matrix formula 



Rather than to compute the limit of the r.h.s. of eq. ( 5.1 6| ) directly, it is easier to compare eq. (1.5) (for some 
trial r-matrix) with eq. ( |6.4[ ) and look at the " discrepancy" (if any) . As a " first approximation" let us try 
the standard classical r-matrix 



2(z-Z- 1 ) & + Z 1 ) / ® / + 2cr l® cr l + 2cr 2®Cr2 + (z + Z : )tT3(80-3] 



(6.5) 



/ z + z- 1 




2 



V 



o \ 




+ Z- 1 J 



which works in the lattice SG model. It is skew-symmetric, i.e. for Z\ ^ z 2 it enjoys the property ( z i/ z 2) 
~ r< 2i { z 2 / z\) where we use the notation from Sect. 5.1. 



Plugging the vectors (5.17) into fll.5| ) with this r-matrix and using eq. (2.21), we compute the "discrep- 
ancy" : 



VVn-A^ (P(l)\r(°\z)\a(l - 1)) 
2A un A pu Ap' M (/3(Z)|a(£ — 1)) 



Mi{kz) 



1 / P'u 











\-l _ \-l 

"pu A un 



(6.6) 



which stands in the right hand side. The result means that our "first approximation" is not so bad. To get 
zero in the r.h.s., we need two corrections - one in the Mi(z) and one in the r^°\ The first one is rather a 
matter of definition. Using the redundant freedom mentioned above, we are free to redefine the M-operator: 

1 



(6.7) 
(6.8) 

■^(^pn^p'n A un) 

The K-term violates the skew-symmetry; nevertheless, the A K%> (z) obeys the classical Yang-Baxter equation 
written in the form modified for not necessarily skew-symmetric r-matrices [p8i : 



Mi(z) = Mt(z) + - \Xp, u + a;J - \ m i ) I . 

The second one is more important. Let r^ K \z) be the modified r-matrix: 

r (K) (z) = r (0) (z) + Kl <g> cr 3 , 



r[2 > (zi/z 2 ) 1 r i £ ) {z 1 /z 3 ) + r^ ) (zi/z 2 ),r^ ) (z 2 /z 3 ) + r$ (z 3 /z 2 ), r[$ (zi/z 3 ) 



= 0. 



(6.10) 



Incorporating the corrections into eq. ( |6.6| ), we arrive at the r-matrix formula 

M t (kz) VVn-AL (P(l)\r (K) (z)\a(l 1)) 



(6.11) 



2A„„A pu Ap'„ (/?(/) | a(/ — 1)) 
At k = (i.e. Ap' U A„„ = Ap„Ap/„) one gets the standard r-matrix. (This just happens for the lattice SG 



model.) Note that in this case £ = 1 + 0(e 2 ) (see (5.2)), so the right classical r-matrix is obtained as a 
" germ" of the quantum one: 

R(z;q,i)^(z^z- 1 )l0l 



r (0) (z) = lim 



o (q - \)[z - z^ 1 ) 
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(provided | = 1 + 0(e 2 )!). 



Remark The non-skew-symmetric r-matrix (6.8) signifies that the L-operator at k ^ is not ultralocal. In 
this case Poisson brackets between elements of the monodromy matrix T are given by a general quadratic 
Poisson bracket algebra 

{7i, T 2 } = r^ 2 TiT 2 + Tisf 2 T 2 - T 2 s^ 2 Ti - T x T 2 r^ 2 

studied (together with its quantum version) in |29|j , |30|] , f3l"l] . The matrices , meet a number of 
consistency conditions which ensure antisymmetry and Jacobi identity for the Poisson bracket. Besides, in 
integrable systems the constraint r + + s + = s" + r~ holds true. Under these conditions, the (non-skew- 
symmetric) r-matrix r + — s~~ satisfies the Yang-Baxter equation of the form ( |6.10| ). In our case the quadruplet 
r + , s + , s~ , r~ is as follows: r + = A°\ s + = —kos ® I, s~ = —nl (8> 03, r~ = + n(I ® 03 — 03 ® I), so 



7 Conclusion 



The main result of this work is the i?-matrix representation ( |5.15 ), (5.16) of the local L-M pair for the 



classical discrete HM model. In our opinion, the very fact that typical quantum i?-matrices naturally arise 
in a purely classical problem is important and interesting by itself. To recover them inside classical discrete 
problems, one has to pass to Hirota's bilinear formalism and make use of Miwa's interpretation of discrete 
flows. As a by-product, we have shown that components of the vectors a), \p) (representing the L-operator 
at the degeneracy point) are r-functions. 

Let us recall that the quantum Yang-Baxter equation already appeared in connection with purely classical 
problems, though in a different context |32|, Q. However, the class of solutions relevant to classical problems 
is most likely very far from R- matrices of known quantum integrable models. In our construction, the role of 
the quantum Yang-Baxter equation remains obscure; instead, the most popular 4x4 trigonometric quantum 
i?-matrix is shown to take part in the zero curvature representation of a classical discrete model. We believe 
that a conceptual explanation of this phenomenon nevertheless relies on the quantum Yang-Baxter equation. 

We should stress that the "quantum deformation parameter" q of the trigonometric R- matrix in our 
context seems to have nothing to do with any kind of quantization. This fact suggests to extend this hidden 
R- matrix structure to the quantum level. In this case the vectors |a), |/3) will have operator components. If 
such an extension does exist, this would mean that there are two different R- matrices in a quantum integrable 
model rather than one. At the same time this would be an intelligent explanation of our .R-matrix formulas 
- in the classical limit one of the i?-matrices degenerates to a classical one while the other one survives. 

As far as other possible projects are concerned, there are many problems apparently arising from this work. 
To mention only few, we point out first of all a detailed comparison with the hamiltonian approach, which 
is missing in our exposition. More specifically, it would be very instructive to identify the Poisson bracket 



algebra for elements of the L-operator (4.12) (in general non-ultralocal) and, starting from this algebra, to 



give an alternative derivation of the r-matrix formula (6.11) for the M-operator. Another interesting question 



is to find a discrete time analogue of the non-local generating function (1.2) of Af-operators. At last, we 
point out the challenging problem to pursue all the program for the " master model" with 2x2 L-Af-pair - a 
discrete analogue of the Landau-Lifshitz model (XYZ magnet). In this case the spectral parameter lives on 
an elliptic curve and the proper extension of the bilinear formalism is not completely clear. (In particular, 
how the reduction condition (|3.2[) looks like?) 
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